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ABSTRACT 

Using the Infrared Spectrograph aboard the Spitzer Space Telescope, we observed multiple epochs 
of 11 actively accreting T Tauri stars in the nearby Taurus-Auriga star forming region. In total, 88 
low-resolution mid- infrared spectra were collected over 1.5 years in Cycles 2 and 3. The results of this 
multi-epoch survey show that the 10 fim. silicate complex in the spectra of two sources — DG Tau 
and XZ Tau — undergoes significant variations with the silicate feature growing both weaker and 
stronger over month- and year-long timcscalcs. Shorter timcscale variations on day- to week-long 
timescales were not detected within the measured flux errors. The time resolution coverage of this 
data set is inadequate for determining if the variations are periodic. Pure emission compositional 
models of the silicate complex in each epoch of the DG Tau and XZ Tau spectra provide poor fits to 
the observed silicate features. These results agree with those of previous groups that attempted to fit 
only single-epoch observations of these sources. Simple two-temperature, two-slab models with similar 
compositions successfully reproduce the observed variations in the silicate features. These models hint 
at a self- absorption origin of the diminution of the silicate complex instead of a compositional change 
in the population of emitting dust grains. We discuss several scenarios for producing such variability 
including disk shadowing, vertical mixing, variations in disk heating, and disk wind events associated 
with accretion outbursts. 

Subject headings: circumstellar matter — infrared: stars — solar system: formation — stars: pre-main 
sequence 



1. INTRODUCTION 

The spatial distribution, size, and composition of sil- 
icate dust grains observed in circumstellar disks sur- 
rounding young pre-main-sequence stars result from 
a variety of physical and chemical evolutionary pro- 
cesses. As they play a key role in the evolu- 
tion of circumstellar disks toward planetary systems, 
an array of observational techniques has been cm- 
ployed to study the characteristics of these dust grains 
(sub-millim eter, millimeter, mid-in f rared, polari z ation , 
etc., (e.g. iWeintraub et all 119891: iStrom et all 119891 : 
iBeckwith et al.1 119901: ISkrutskie et all I1990D . Building 
upon much ground-based infrared photometry and imag- 
ing sensitive to the blackbody emission and scattered 
light from dust grains in circumstellar environments, the 
Spitzer Infrared Spectrograph (IRS) has amassed an ex- 
tensive catalog of near- to mid-infrared spectral energy 
distributions (S EDs) of young sources in the nearest star 
forming regions (IKessler-Silacci et al.ll2005l; iFurlan et all 
120061 ISargent et all 120091 : iWatson et al J 120091 ) . Model- 
ing of the SEDs of this large sample of T Tauri stars 
(TTS) — known to possess thick and actively accreting 
disks — succeeds in defining the spatial and size dis- 
tributions of the dust grains in these systems, poten- 
tially giving astronomers direct insight into the evolu- 
tionary status of these planetary building blocks. For 
instance, SED models suggest the presence of 'puffed- 
up' inner disk rims near the dust sublimation radius for 



some sources (INatta et all 120011: IDullemond et al.l 120011: 
IMuzerolle et all I2003L 12001 " and constrain the sizes of 
the optically thin inner disk holes of t ransitional objects 
finding some to be as large as 40 AU (|Strom et al.lll989l: 
I Skrutskie et all 119901 : iCalvet et aLll2002t lEspaillat et al l 
l2007allbf ). 

In addition to the information encoded in the SEDs, 
grain sizes and composition have been inferred from 
mid-infrared emission features formed by the super- 
heated dust grains located at small radii and in the up- 
per atmospheres of these protoplanetary disks. High- 
sensitivity, low-resolution spectral surveys of the 10 fxm 
silicate complex have provided the opportunity for in- 
depth analyses of compositional variations of the dust 
grains with the most detailed models including labo- 
ratory measur ed opacities for a variety of sizes and 
species (see e.g.lBouwman et all2 003: Mc cus et al.ll2003l: 



Honda et all 120061 : iSchegerer et all 120061 : ISargent et all 
20091) . The parameters returned by these models in- 
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elude dust temperatures and mass fractions for differ 
ent grain types, sizes, and structures (i.e., crystalline 
and amorphous, small and large, forsterite, enstatite, 
olivine). Less detailed compositional analyses have been 
conducted by binning the 10 /zm silicate complex into 
various wavelength regions corresponding to the peak 
emission of specific dust grain species. The integrated 
fluxes over these wavelength bins are used to com- 
pute ratios commonly referred to as crystalline indices 
that give an approximation of t he relative abundances 
ivirgijiffiiede g-. Bouwman et al.l 120011: iPrzygodda et al.l 2003 : 



van Boekel et all 120031: IKessler-Silacci et al.l 120051 12006 



Watson et al.H2009ft . 



Several groups have published single-epoch Spitzer IRS 
surveys of TTS in which they note source-to-source mor- 
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pholo g ical differences of the si licate f eature (|Forrest et al.l 
20041: iKessler-Silacti et~aTl 120051; iFurlan et all 120061; 
Sargent et al.l 120091 ; I Watson et al I l2009f ). iForrest et all 
(|2004f l were the first to note a relationship between 
the crystallinity versus amorphous nature of the silicate 
emission feature with the evolutionary status of the dust 
grain population and the disk. The presence of crys- 
talline grains implies that the once amorphous interstel- 
lar grains have been annealed and cooled, consistent with 
heating i n the inner regions of an ev olvi ng protoplane- 
tary d isk. [Kcsslcr- Silacci et alJ (|2005l ) and I Sargent et all 
(|2009f ) concluded the same, suggesting that a disk evo- 
lutionary sequence can explain the variety of shapes ob- 
served for the silicate features in their Short-Low (5.2- 
14.5 urn) and Long-Low (14.0-38.0 fim) IRS spectra of 
TTS. 

The conclusions about the evolutionary status of the 
disks are based solely on single-epoch observations and 
do not account for the dynamic nature of the TTS 
sources, which are characterized by multi-wavelength 
variability on timescales as short as days. Variations of 
the morphology of the silicate emission feature on the 
timescales of days, months, or even years would likely 
complicate such compositional studies of dust grain pop- 
ulations such as those leading to the development of a 
disk evolutionary sequence. 

In order to test the validity of the conclusions drawn 
from the single-epoch studies, we used Spitzer IRS to 
conduct a spectroscopic variability study of the 10 fim. 
silicate emission feature in 11 actively accreting classical 
TTS found in the Taurus- Auriga star forming region. In 
2 of the 11 sources, DG Tau and XZ Tau, significant vari- 
ability of the 10 fj,m feature is observed. The unexpected 
short timescale variations and potential scenarios for this 
behavior are presented in this Letter. A long paper pre- 
senting more detailed modeling of these features and the 
spectra of the sources with no detectable variability is 
forthcoming (Leisenring et al., in preparation). 

2. OBSERVATIONS AND DATA PROCESSING 

Low resolution (R ~ 100) spectra of our selected sam- 
ple of actively accreting TTS were obtained using IRS 
in Short-Low mode (3'.'6 slit and 5.2-14.5 /im wave- 
length coverage) in Spitzer Cycles 2 and 3. In order to 
probe silicate variability on week-, month-, and year- long 
timescales, each source was observed twice in each of the 
two approximately month-long visibility windows avail- 
able per Spitzer Cycle. As a result, we collected a total of 
88 spectra spanning a baseline of 1.5 years 3 . Archival IRS 
observations have been included to extend the temporal 
coverage to a longer timescale. Table [T] summarizes the 
multiple epochs of observations of DG Tau and XZ Tau 
and includes measurements for the continuum color and 
the integrated flux of the silicate complex. 

The spectral images were reduced using the SMART 
software package cr eated by the IRS instrument team 
(iHigdon et al.ll2004D . the Basic Calibrated Data (BCD) 
produced by the IRS data pipeline, and routines stan- 
dardized b y the Guaranteed T ime Observer IRS Disks 
team (e.g.. ISargent et aTll2006f ). Due to the broad wave- 

3 For a couple of sources, only one observation was made in the 
first visibility window of Cycle 2 and three were made in the other 
for a given cycle. 
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Fig. 1. — Examples of the wavelength-dependent correction fac- 
tors employed to account for flux loss due to the target being offset 
from the center of the slit by 0'.'5 (solid line) and l'.'O (dashed line). 
The blue portion and red portion of the spectra corresponds to 
Short-Low order two (5.2-7.7 fim) and one (7.4-14.5 /im), respec- 
tively. The green portions of the spectra correspond to the bonus 
order segment (7.3-8.7 fim). Error bars are based on a pointing 
accuracy of 0'.'2 (SSC Helpdcsk, private communication) and rep- 
resent the standard deviation in the absolute flux corrections to 
the spectra and not errors in flux from one wavelength resolution 
element to the next. 

length coverage of IRS, the shapes and fluxes of the con- 
tinuum and solid-state features can be affected by vari- 
ations of the pointing of the telescope for different nod 
positions and from epoch to epoch. Using an average 
of well-pointed staring observations (AORs: 16344064, 
19438080, 20957184, 22165760, 24240896, 24574464^ 
27392768, and 28709632) and spectral maps (AORs: 
16295168 and 19324160) of the standard star HR 7341, 
the slit offsets were reconstructed and the science spectra 
were corrected for any epoch-to-epoch mispointings 4 . Af- 
ter applying the correction factors corresponding to the 
appropriate mean slit offsets (see Figure [Hand Table 1, 
Column 6), we confidently rule out slit offsets as a source 
for the variations observed in the shapes and fluxes for 
the silicate features reported in Section 3. In addition, 
given the size of the slit (3'.'6 for IRS), slit offsets, and 
the T/2 to 3'.'6 FWHM of the point spread function for 
Spitzer over the relevant wavelength range of 5 to 15 /im, 
the spectral variability o bserved in the XZ Tau binary 
system (separation of 0'.'3: lHaas et al.l (|1990f )) is unlikely 
to be caused by light loss from either of the stellar com- 
ponents. 

The spectra with continuum fits in Figure [2] illustrate 
the technique employed for removing the continuua from 
the target star spectra and demonstrate that the ob- 
served variability of the silicate feature is independent of 
variations in the long and short wavelength continuum 
regions. In order to quantify the shape of the contin- 
uua and underlying dust continuum emission, we defined 
two wavelength bands centered at 6.0 and 13.5 fim with 
widths of 0.8 and 1.0 /im, respectively, and the color in- 
dex ([6.0] - [13.5]; see Table 1). 

3. RESULTS 

4 The correction factors applied to each spectrum, along with the 
details of the reconstructed pointing measurements will be included 
in a Leisenring et al. (in preparation). 
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TABLE 1 

Observations, Continuum and 10 fxM Feature Measurements 



Star 


AOR 


Date Observed 


SL1 


& 


SL2 


Mean Slit 


[6.0] - [13.5] 


Silicate Flux 


Name 


ID 


Gregorian 


Julian 


Cycles 


X 


Ramps 


Offset (") 


Color 3 


(10- 14 W m~ 2 ) b 


DG Tau 


3530496 


2004-03-02 


2453067 


2 


X 


6 sec 


-1.15 


2.63 ± 0.007 


1.93 ± 0.370 




15110144 


2005-10-09 


2453653 


8 


x 6 sec 


-1.97 


2.80 ± 0.009 


4.05 ± 0.908 




15115008 


2005-10-15 


2453659 








-2.15 


2.87 ± 0.009 


3.99 ± 0.866 




15115264 


2006-03-17 


2453812 








0.00 


2.90 ± 0.011 


6.24 ± 0.076 




15115520 


2006-03-22 


2453817 








0.00 


2.90 ± 0.011 


6.00 ± 0.077 




19488000 


2006-10-17 


2454026 








-0.23 


2.66 ± 0.010 


5.30 ± 0.057 




19487744 


2006-10-22 


2454031 








-0.40 


2.64 ± 0.008 


3.98 ± 0.166 




19487488 


2007-03-21 


2454181 








0.23 


2.75 ± 0.010 


5.78 ± 0.102 




19487232 


2007-03-28 


2454188 








0.96 


2.85 ± 0.011 


6.52 ± 0.516 


XZ Tau 


3531776 


2004-03-04 


2453069 


2 


X 


6 sec 


-0.30 


2.35 ± 0.008 


4.43 ± 0.163 




15111680 


2005-09-12 


2453626 


8 


X 


6 sec 


0.46 


2.37 ± 0.006 


1.29 ± 0.087 




15118336 


2006-03-09 


2453804 








0.36 


2.42 ± 0.007 


4.16 ± 0.090 




15118592 


2006-03-17 


2453812 








0.66 


2.44 ± 0.009 


4.04 ± 0.212 




16875776 


2006-03-21 


2453816 








0.75 


2.42 ± 0.006 


3.67 ± 0.238 




19494144 


2006-10-17 


2454026 








-0.32 


2.50 ± 0.007 


6.12 ± 0.094 




19493888 


2006-10-22 


2454031 








-0.57 


2.50 ± 0.007 


5.86 ± 0.236 




19493632 


2007-03-22 


2454182 








0.26 


2.34 ± 0.010 


3.92 ± 0.066 




19493376 


2007-03-28 


2454188 








0.10 


2.34 ± 0.008 


3.97 ± 0.046 



Errors for the color temperature are based off of relative flux uncertainties. 

Errors for the flux of the 10 /mi silicate feature are based off of pointing-induced flux uncertainties. 
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Fig. 2. — Two epochs of spectra for both DG Tau and XZ Tau are 
the solid black lines in each panel. The smooth gray line represents 
the best third-order polynomial fit to the short (5.2—6.4 fim) and 
long (13—14 fim) wavelength portions of the spectra. For compari- 
son, the residual spectra (black dot-dashed lines) depict the shapes 
and relative fluxes of the silicate features, clearly illustrating the 
large variability observed between these epochs of observations. 

Of the 11 sources included in our program, only 
DG Tau and XZ Tau exhibit variability of the sili- 
cate complex over month- and year-long timescales. 5 In 
both systems, the silicate emission feature is observed 
to weaken and strengthen over the three years worth 
of observations, with the short wavelength portion near 
9 ^m in XZ Tau falling nearly to the level of the con- 
tinuum. The silicate variability in th e spectrum of 
DG Tau has previously been observed (IWooden et all 
l2000t IWoodward eFall l200i ISitko et al l 120081) and is 
confirmed by these observations, while the variations in 

5 We note that the preliminary data analysis found evi- 
dence for short day-t o-day and week-long timescale variations 
(Lciscnring et al. 2007). However, once the correction factors for 
mispointings were applied, these short timescale variations were no 
longer confidently detected. 



the spectra of XZ Tau are observed for the first time. 

Using the JD 2 453069 observation of XZ Tau, 
ISargent et al.l (|2009f ) determined that the silicate fea- 
ture is dominated by amorphous grains by fitting a 
two-temperature blackbody with solid angle dependence 
and laboratory measured opacities for a variety of dust 
grain sizes and compositions. By inspection of the 
XZ Tau spectra (Figure [3]), for all epochs, this source 
lacks evidence of emission from crystalline forsterite 
at 11.3 /im, a common signature of sources deter- 
mined to poss ess relatively high abundances of crys- 
talline grains (.Forrest et all 1200 4; Kcsslcr- Silacci et al.l 
120061: ISargent et al.ll2009D . On the other hand, the sil- 
icate emission feature for DG Tau appears dominated 
by crystalline forsterite in all epochs of observations pre- 
sented here 6 . Detailed modeling of the silicate feature 
is required to extract the crystalline-to-amorphous ratio 
and will be presented in Leisenring et al. (in preparation) . 

In Figure 21 the silicate flux appears to be marginally 
correlated with the [6.0] — [13.5] color index for both DG 
Tau and XZ Tau. However, such weak correlations do 
not provide compelling evidence for a common physical 
mechanism that is directly responsible for the variations 
in the silicate feature and also for changes in continuum 
emission. We hope to further inquire into this relation 
with (near-)simultaneous multi-wavclcngth observations. 

3.1. Preliminary Modeling 

In addressing the evolutionary sequence previously 
proposed as an explanation for the amorphous and 
crystalline silicate abundances, we adopted the two- 
temperature dust models (inclu ding the sam e dust 
species and related opacities) of ISargent et al.l (]2009f ) 
(see Equation 1) and performed compositional fits to the 
10 /urn silicate feature at each phase of modulation for 
DG Tau and XZ Tau. The models returned statistically 
significant fits only for the epochs possessing maximum 

6 IWatson etaLl |(2009T ) were unable to confidently extract a 
meaningful crystalline index from a single-epoch observation for 
this source. 
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Fig. 3. — All epochs of spectra acquired for both DG Tau and XZ Tau with Spitzer IRS. In each case, the silicate feature represents 
a different compositional signature. The silicate feature for DG Tau appears strongly peaked near 11.3 (im the location of crystalline 
forsterite emission, while the very smooth silicate feature of XZ Tau peaks at shorter wavelengths and lacks this notable crystalline 
emission. Therefore, by inspection and the previously used crystalline indices, the silicate dust grain populations are characterized by an 
abundance of crystalline dust grains in the disk of DG Tau and amorphous grains for XZ Tau. Note that the absolute flux error bars show 
that the small variations are not likely to be real, which rule out week-long silicate variability in these observations. 

silicate feature do not necessarily indicate a change in 
the composition of the emitting grains. The introduc- 
tion of an absorbing column of dust grains explains the 
inability of previous silicate decomposition models to fit 
silicate emission features in DG Tau and other similar 
sources (jSargent et alJl200ft IWatson et al.ll2009T >. 

4. INTERVENING COOLER DUST GRAIN SCENARIOS 

The timescales over which the variations occur strongly 
suggest that the phenomena associated with intervening 
cooler dust must occur on dynamical timescales consis- 
tent with the motion of dust in the disk at radii less 
than or equal to 1 AU. These short timescales rule out 
a clumpy dust envelope as the source for the interven- 
ing column of cooler dust. Th i s infer ence agrees with 
the findings of Woodw ard et al.l (|2004f ). in which a one- 
dimcnsional radiative transfer model showed that an 
envelope capable of producing the variability they ob- 
served in DG Tau would also completely obscure the star. 
Therefore, we believe that it is unlikely that the inter- 
vening cooler grains are associated with remnants of an 
envelope. 

In an attempt to explai n day- to week-long si licate 
and continuum variability, Muzcrol le et alj (|2009t ) dis- 
cuss a scenario of disk shadowing in which a warped, 
puffed up region of the inner disk rim casts a shadow 
across a portion of the disk surface. There is potential 
for this phenomena to explain the longer timescale vari- 
ability presented here (|Sitko et al.l [2008). In addition, 
the recent discovery of a third companion located within 
0'/09 (13 AU) of the southern component of XZ Tau 
(jCarrasco-Gonzalez et al.ll2009T l increases the possibility 
for variations in disk illuminatio n as previously described 
by IWatson fc St apclfcldtJ (|2007l ) in reference to the vari- 
ability they observe in the HH 30 system. 

In additio n to disk shadowing and illumination, tu rbu- 
lcnt mixing ([Turner et al.ll200a [Balsara et alJl2008h and 
the presence of dusty disk winds may also provide rea- 
sonable explanations of the observed variability. Turbu- 
lent mixing in the disk capable of dredging cooler dust 
up to the surface of the disk may provide a plausible 
mechanism for varying the strength and nature (emis- 
sion versus self-absorption) of the silicate feature. Altcr- 



[6.0]-[13.5] 

Fig. 4. — Silicate flux plotted vs. the [6.0] — [13.5] color for both 
DG Tau and XZ Tau. We obtain respective Pearson correlation 
coefficients between these two parameters of 0.61 and 0.66 for 
DG Tau and XZ Tau. The probabilities that these correlations 
could h ave been pro duced from a random distribution (see, for ex- 
ample, Taylor 1997) arc 8.1% and 5.3%, respectively. Thus, we 
find only a weak correlation between the strength of the silicate 
feature and the continuum flux. 

silicate emission (DG Tau: JD 2453653 and 2453659; XZ 
Tau: JD 2454026 and 2454031). Therefore, the models 
were unable to provide a reliable fit to the silicate feature 
when the emission was diminished or self-absorbed. The 
failure of this simple method led us to develop a two- 
temperature two-slab model of the dust emission from 
the surface of the disk. In order to produce the ob- 
served variations of the silicate feature, the top, interven- 
ing layer was assumed to be cooler than that of the lower, 
emitting region. This two-slab model confidently fits all 
epochs of the silicate emission features for both sources, 
especially when both slabs possessed dust grains with 
similar compositions. The results of this model 7 sug- 
gests that the observed variability can be explained by 
an intervening population of similarly composed cooler 
dust grains. Therefore, the variations observed in the 

7 The models and results will be detailed in Leisenring et al. (in 
preparation). 
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na tively, the presence of dusty d isk winds as described 
bv lTambovtseva fc Grinml (|2008l ) may be capable of loft- 
ing material above the disk surface providing a cooler 
veil of similarly composed dust grains through which to 
view the warmer, emitting dust; hence, producing self- 
absorption of the emitting silicate region. In both sce- 
narios, the vertical mixing and disk winds must be vari- 
able to produce the inhomogeneities required to explain 
the observed variability of the silicate feature. Given 
the variable nature of the accretion processes in these 
sources, fluct uating disk winds and vertical mixing is to 
be expected (|Barv et al.ll2008HNguven et al.ll2009f ). 

5. CONCLUSIONS 

The 10 /im silicate emission complex observed in the 
mid-infrared spectra of classical TTS and formed by 
super-heated dust grains in the upper atmospheres of 
the innermost regions of their disks can vary on month- 
long timescalcs. Better temporal resolution is required to 
determine if the variability is periodic. Our spectral de- 
composition models were incapable of fitting the silicate 
feature at every phase of modulation, suggesting that 
variability may be an explanation for some of the diffi- 
culties previous groups have encountered with modeling 
single-epoch spectra of individual stars. The simple two- 
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slab two-temperature model suggests an optical depth 
effect may explain the observed variations. Since the in- 
tervening cooler material may be of similar composition, 
the variability may not directly impact the disk evolu- 
tionary sequence previously constructed for TTS, based 
on amorphous versus crystalline abundances and dust 
grain sizes. However, we note that these models require 
a number of parameters and the question of degenerate 
solutions needs to be addressed. Additional multi-epoch 
observations of these sources and others are required to 
determine the fraction of disks with variable silicate emis- 
sion and provide more variations to model, helping to 
constrain acceptable parameters. With the loss of IRS on 
Spitzer, this project will turn to ground-based observa- 
tions with less sensitivity, but still the ability to measure 
the previously observed silicate variability. 
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